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O '. ABSTRACT 
(N 

{^JqI Context. Herbig Ae/Be stars are intermediate-mass pre-main sequence (PMS) stars showing signs of intense activity and strong 

, stellar winds, whose origin is not yet understood in the frame of current theoretical models of stellar evolution for young stars. The 

evolutionary tracks of the earlier Herbig Ae stars cross for a significant fraction of their evolution towards the main sequence the 
theoretical PMS instability strip located roughly in the same area of the HR diagram as the 6 Scuti variables. Many of these stars 
^■f^ ' exhibit pulsations of 6 Scuti type. 

Aims. We carry out a thorough analysis of the line profile variations of the pulsating prototype Herbig Ae star HD 104237 based 
on high-resolution spectroscopical time-series, identifiy the dominant non-radial pulsation mode, and perform the most accurate 
determination of its fundamental parameters in preparation of an asteroseismic modeling. 
' Methods. HD 104237 is a pulsating Herbig Ae star with eight detected frequencies based on the analysis of radial velocity variations. 

. , In this article, we reinvestigated an extensive high-resolution quasi-continuous spectroscopic data set in order to search for very 

. faint indications of non-radial pulsations in the line profile. To do this, we worked on dynamical spectra of equivalent photospheric 

Qh' (LSD) profiles of HD 104237. A 2D Fourier analysis (F2D) was performed of the entire profile and enabled the identification of the 

dominant non-radial pulsating mode. The temporal variation of the central depth of the line was studied with the time-series analysis 
tools Period04 and SigSpec. The development of a highly accurate continuum normalization method enabled us to determine reliably 
the equivalent width of chosen absorption lines and to perform a new determination of the fundamental stellar parameters. 
Results. Following the previous studies on this star, our analysis of the dynamical spectrum of recentered LSD profiles corresponding 
to the 22"'' - 25"" of April 1999 nights spectra has confirmed the presence of multiple oscillation modes of low-degree C in HD 104237 
and led to the first direct detection of a non-radial pulsation mode in this star: the dominant mode Fl was identified by the Fourier 2D 
method having a degree t value comprised between 1 and 2, the symmetry of the pattern variation indicating an azimuthal order of + 
^ ' 1. The detailed study of the fundamental stellar parameters has provided a T^a, log g and iron abundance of 8550 + 150 K, 3.9 ± 0.3 

■ and -4.38 ± 0.19 (i.e. [Fe/H]= -1-0.16 ± 0.19 ), respectively. 
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1. Introduction lines are observed and P Cygni profiles of Ha, Hj3, Mg II h 

and k lines present in the spectra of some of these stars indi- 

Pre-main sequence (PMS) stars of intermediate mass (2 to 10 ^^^^ j^e presence of strong and structured stellai- winds (see e.g. 

Mo) are known as Herbig Ae/Be (hereafter HAeBe) stars, fol- praderie et al. 1982, 1986; Catala & Talavera 1984; Catala et d] 

lowmg a first classification by Herbig (1960). As their name indi- 19363 b; Catala & Kunasz 1987; Catala 1988t lCatala et al.ll99ir 

cates, they ai-e of spectral type A or B and luminosity classes III Bohm & Catala.J995; Bohm et al 1996) . 

to V (for a review, see Wate rs & Waelkenslll998l) . Their Spectral 

Energy Distribution (SED) is characterized by infrared (IR) ex- a wide-spread idea consists in invoking magnetism as being 

cess and ultraviolet (UV) extinction mainly due to thermal re- responsible for such active stellar phenomena, requiring either 

emission from UV to IR of a circumstellar (CS) dust shell or presence of a primordial magnetic field or a dynamo mecha- 

disk (or both), making the photometric assessment of spectral nism whose basic ingredients are classically convection and ro- 

types and stellar temperatures of these stai's rather uncertain. ^^i^qj^ Qn one hand, recent specti-opolarimetric observations of 

HAeBe stars show signs of intense stellar activity, variabil- Herbig stars indicate the presence of significant magnetic fields 

ity and winds. Indeed, the frequently observed emission in Mg only in a small 10% fraction of them (Wade et al. 2007), these re- 

II h and k lines, Ca II IR triplet, Ca II H and K line. He I suits being in agreement with the primordial fossil field hypoth- 

5876 A, many Fell lines, Na I D, N V, Si IV and C IV reso- esis. On the other hand, considering a possible dynamo mecha- 

nance lines, imply the high temperature of a chromosphere or nism, the position of HAeBe stai-s in the HR diagram indicates 

corona to be formed. Short-term vai-iability of many of these that they ai-e in the radiative phase of their quasi-static contrac- 

tion towards the main sequence (Iben 1965; Gilliland 198^. 

* Based on observations collected at the 1.9m SAAO Radclifi'e tele- Therefore, in absence of subphotospheric convection, a classi- 

scope cal solar-type magnetic dynamo mechanism can not be at work. 
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The HAeBe stars activity remains therefore quite paradox- 
ical in the frame of current theoretical evolutionary models for 
PMS stars, although there exist so me approaches to understand 
its origin. As possible explanations. iPalla & Stahleil(ll99 0) sug- 
gested that out er an convection zone could be induced by deu- 
terium burning. [Vigneron et alj (l990) proposed the internal ro- 
tation as energy source, the friction exerted at the stellar surface 
by the angular momentum losses creating a non-convective tur- 
bulent layer, r esponsible for a dynam o effect; this model was 
later refined bv lLignieres et al.l (Il996l) . 

An external origin, involving an accretion disk dissipating 
gravitational energy in to a boundary layer at the stellar surface, 
was invoked bv Bertout et al.l (Il988l) to explain activity in classi- 
cal TTauri stars, TTauri star s being the lower mass P MS counter- 
parts to HAeBes. However, iBohm & Catalal ( Il995h showed that 
the energy flux emitted in winds and some activity tracers seems 
to increase with effective temperature, which is in favor of an 
internal activity source rather than a CS origin. Thus, as of to- 
day, growing evidences tend to indicate that the energy needed 
to produce this activity might be of internal stellar origin, but 
no definite answer has been provided. It is a major concern for 
testing young stellar evolutionary theory to solve this still open 
question about HAeBe stellar activity, by constraining the inter- 
nal structure of these objects using asteroseismic techniques, i.e. 
the analysis and modeling of stellar pulsations, if observed. 

Photometric and spectroscopic variability have indeed been 
observed in HAeBe stars, on very different time scales (from 
minutes to years). Large and irregular photometric variations 
(> O^.S) with time scales from weeks to years are thought to 
be due to variable obscuration by CS pat chy dust clouds, par- 
ticula rly strongly in the younger objects (Ivan den Ancker et al.l 
Il998h . Variations of the order of a tenth of magnitude (mag) in 
amplitude, with time scal es of hours to days, can be explained as 
due to clumpy accretion (Ivan den Ancker elTal . 1998 and refer- 
ences therein) or due to chromospheric activity, including winds, 
modulated by the star's rotation ( Catala et al..,1999.) . Stellar pul- 
sations, which are here the subject of our interest, induce typ- 
ically a millimagnitude (mmag) photometrical variability with 
time scales of minutes to hours. 

The existenc e of su c h stel lar pulsations in HAeBe stars 
are known since i Bregerl (Il972h discovered two pulsating can- 
didates, V588 Mon and V589 Mon, in the young open clus- 
ter NGC 2264. Subsequently, two additional pulsating Herbig 
stars were detected: H R 599 9 (Kurtz & Marang 1995|) and 
HP 104237 jP onati et aLl ll997h . Motivated by these detections, 
iMarconi & Palla (J998i) investigated theoretically a PMS insta- 
bility strip for the first three radial modes of pulsations and 
concluded on a topology of this instability strip, corresponding 
roughly to the same area in the HR diagram as the 5 Scuti vari- 
ables. The study of iMarconi & PaUal (Il998h showed that most 
of the Herbig Ae stars are expected to cross it for a significant 
fraction of their evolution to the main sequence (5 to 10% of 
their PMS phase). Since the precise location and boundaries of 
the PSM instability strip has not yet been observationnally con- 
strained, it is necessary to identify and study the largest number 
of PMS pulsating objects, giving priority to multiperidodic pul- 
sators which are suitable and promising candidates for future as- 
teroseismic modeling, and, by this mean, supplying constraints 
on theoretical evolutionary pre-main sequence models. 

As of today, a significant number of Herbig Ae stars, either 
field stars or members of young open clusters, have indeed re- 
vealed to be pulsating at timescale typical of 6 Scuti stars, i.e. 
with short periods (from a; 20 min to several hours) and small 
amplitudes (from mmag to few hundreths of mag in case of pho- 



tometry, and less than typically 1 - 2 kms for radial veloc- 
ity studies) (for reviews see e.g. 'Catala"2003t IMarconi & PaUal 
2004; Zwintz et al. 2004; Zwintz 2008). Most of them have been 
stud ied in photom etry, but only very few in spectroscopy (e.g. 
Bohm et al. I I2OOI 12009 ). However, due to cancellation effects, 
only modes of low degree are detected in photometry or from 
radial velocity curves. Therefore, high-resolution spectroscopy 
is necessary to study the weak variations induced by non-radial 
pulsations in the rotationally-broadened photospheric line pro- 
files, and thus to carry out a comprehensive pulsational analy- 
sis in preparation for seismic modeling. It should be noted that 
binarity is commo n amongst HAeBe stars with approximately 
60% (Baines et al.' 2006). As we will see in Sect.|2 HD 104237 
is also a spectroscopic binary, which complicates tremendously 
the analysis of its spectrum. 

This article is the second of a series on the prot otype Herbig 
Ae star HD 104237, starting with the article by iBohm et aU 
(2004). The first article presented the detection of multiperiodic 
oscillations in the radial velocity curves of the star and the anal- 
ysis of its binary orbit, amongst others. In this present article, 
we describe the direct search and detection of oscillations in the 
line profiles, the mode identification of the dominant non-radial 
pulsation mode, as well as a fundamental stellar parameter de- 
termination by spectroscopic means. This work represents there- 
fore the basis for the asteroseismic mod eling we will present in 
a forthcoming article dFumel et al.ll201 ll) . 

The article is structured as follows: in Sect. |2] we review in 
detail previous works on HD 104237, Sect.|3]sums up the spec- 
troscopic observations we worked with and data reduction. Sect. 
|4| presents our detection and identification of non-radial pulsa- 
tions in HD 104237, Sect. |5] describes the spectral continuum 
determination tool we have developed to obtain optimized high- 
quality spectra, in Sect. |6] we detail the fundamental parameter 
determination of HD 104237 we carried out and finally we con- 
clude and discuss our results in Sect. |7] 



2. Previous works on HD 104237 

IHu et al.l (Il989l Il99ll) derived from photometric and spectro- 
scopic studies that HD 104237, or DX Cha, is a member of 
the Herbig Ae/Be group. Moreover, it is a prototype of this 
group of stars and a particularly suitable target since it is very 
bright (njv = 6.6), which enables to observe it in high reso- 
lution spectroscopic mode. It is located in the e Chameleontis 
young stellar group at a H ipparcos distance of 116 + 8 pc 
(Ivan den Ancker et al.1 1 1997b . HD 104237 is actually a multi- 
ple system whose primary component (hereafter simply called 
HD 104237 or the "primary") is a Herbig Ae star with sev- 
eral low-mass compan ions at separations between 1" and 15" 
(iFeigelson et al. I l200l . including a close K3 companion in an 
eccentric 19.9-day orbit (hereafter HD 104237b or the "sec- 
ondary"). This very close component of approximate spec- 
tral class K3 forms with the primary a spectroscopic binary 
(iBohm et al.l 12004) . Indeed, K3-type spectral features (includ- 
ing Li I 6707 A and Ca I 6718 A li nes) are observed in the 
primary spectra ("Feigel son et al.l l2003). that makes its spectro- 
scopic analysis trickier since one has to take into account the pol- 
lution generated by the faint secondary component in the domi- 
nant primary spectrum. 
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2.1. Previous fundamental parameter determination 

The strong IR excesses and the UV-extinction due to the pres- 
ence of CS dust , the emission components in many absorption 
Unes, the P Cygni profiles observed notably in the Balmer lines 
and the spectroscopic binary pollution make the determination 
of most fundamental stellar parameters of HD 104237 difficult. 
Thus, the spectral types avalaible in the literature for HD 104237 
range from AO to A8 (implying a large range in possible T^ff 
values) and its bolometric luminosity from » 20 to 60 Lq. The 
associated stellar parameters such as the mass and the age of the 
star can be estimated from temperature and luminosity with stel- 
lar evolutionary models. They rely therefore strongly both on 
the accuracy of the location of the star in the HR diagram and 
on the adopted evolutionary models and are also not precisely 
determined, as well as all the parameters that result from them 
(e.g. log g). The most important parameters of HD 104237 (pri- 
mary and secondary) found in the literature are chronologically 
summarized in Table [1] and are detailed below. 

iHu et alJ ([1989 ) based their first determination of a AOVpe 
spectral type on the absorption profiles of the Balmer lines and 
the strong Ca II K line (3933.7 A). Subsequently, the presence 
of emission components in the Ha and H/S was understood to be 
mislea ding the determ ination of the spectral type. In a second 
article EuetaDdMIl) 

therefore derived an A4Ve spectral type 
from UV absorption lines. From an integrated radiant flux (0.1 
to 10 fim) of 8.62 x 10"** erg cm"^ s"' and a distance of 88 pc 
they inferred a log (L^/Lg) of 1.3 4. Using the dwarf tempera- 
ture scale of Schmidt-Kaler, (|1982|) . they converted this spectral 
type to a of 8450 K from whi ch they obtained a stellar ra- 
dius of 2.2 R7^. lBrown et alJ (Il997h carried out a study of the hot 
disk wind of HD 104237 by means of UV spectra and assumed it 
to be of spectral type A7IVe. Using the accurate absolute astro- 
metric and photometr i c data yielded by the Hipparcos satellite, 
Ivan den Ancker et al.l (11997') computed new Tetj, log (L^/L©), 
mass and age values for HD 104237 (cf. Table[TJ. Although their 
method included corrections taking into account the fact that 
HAe/Be stars often exhibit anomalous extinction law for their 
circumstellar material and UV spectra often lead to overestima- 
tion of temperatures due to the presence of heated layers in the 
immediate surrounding of the star's p h otosphere, they u sed as 
a spectral type the non-corrected result o f lHuetalJ(ll989h based 
on optical spectra, whose many lines are polluted by emission 
components. In a next article van den Ancker et al. ( 1998) pro- 
posed a cooler spectral type (A4IVe), in agreement with lHu et al.l 
( 1199 lb . and derived corresponding new parameters. The uncer- 
tainties on the effective temperature were estimated values and 
not determined in a statistical procedure. 

In their invest igation of the environment of HD 104237, 
iGradv et aP (l2004l) derived a new spectral type of A7.5-A8Ve 
based on the comparison between FUV and UV spectra of 
HD 104237 and others stars with well-known spectral types. 
Again, the error bars on the spectral type were empirically de- 
te rmined. This cooler spectral type, close to the determination 
of iBrown et al.1 (Il997h . and the log (L*//.©) derived by the au- 
thors led to an age of 5 Myr, which is noticeably greater than 
the previou s age estimations (van den Ancker et al. 1997, 1998; 
iFeigelson etal.il2003h but cl oser to the ages of flie HD 104237 
quintet low-mass members dFeigelson et al Jl2003l). Usin g evo- 
lutionary tracks and isochrones from IPalla & Stable r (I200I 
and b oth t he stellar parameters from van den Ancker et alJ 
(119981) and iGradv etil. (2004). Bohmetal. (2004J concluded 
on an age of 2 Myr for HD 1 042 37, in agreement wit h 
Ivan den Ancker et al ] d 199711 19981) and iFeigelson et all (l2003l) . 



and independently concluded on a similar value of the stellar 
mass Mi, - 2.2 + O.IMq as previous studies. Moreover they cal- 
culated a mass ratio Mp/Ms of 1.29 + 0.02, based on measures 
of the primary and secondary radial velocities with respect to 
the systemic velocity of the binary system, that enabled them 
to estimate the mass of the secondary to be of 1.7 + O.IM©. 
Assuming the same age for both components, they confirmed 
a K3 spectral type for the secondary and assessed a temperature 
of about 4750 K and a luminosity of one tenth of Lp . Starting 



from t he spectral type proposed bv'Gradv et al.l (|2()04l) 



LuhmanI 



(|2004|) computed a new T^g (with Schmidt-Kalei 11982 temper- 
ature scale) and log {Li^lLo) (from /-band magnitudes). Placing 
these new parameters in a H R diagr am with evolutio nary tracks 
from lPalla & Stahled (1 19991) enabled jLuhmai ([20041) to confirm 
an age of about 2 Myr. Finally, based on low-resolution spectra 
and going by the fact that the strengths of Ca I I and H^ lines un- 
dergo a reverse relationship with temperature, iLvo et al] (l2008l) 
managed to constrain the spectral type of HD 104237, even in 
case of contamination by emission of these lines. Since they ob- 
tain a spectral type of A4 from the Ca II line and A5 from the 
H(5 line, they adopted a spectral type of A4/A5 for this star. 



Acke & Wae 



kens 



(12004 used the T^f( and log g presented 
in Mee us et al.l (12001!) (and references therein) to carry out a 



chemical analysis of several PMS stars including HD 104237. 
They measured the equivalent widths of many chemical ele- 
ments a nd converted them to abundances using the program 
MOOG (ISnedenlll973"l) . but taking into account neither the bi- 
narity of this star (which adds a continuum veiling along with 
additional lines from the secondary) nor the asymmetries or the 
emission components in the line profiles. Using the solar abun- 
dances by lAnders & Grevessd (Il989i) they found a solar-like 
global metallicity ([M/H]= H-0.06 + 0.05) and a solar abundances 
for several chemical elements including Fe ([Fe/H]=-i-0.09+0.19 
where [Fe/H]=log(^) -log(^)). 

L ' J °V A'h ^HD 104237 = V Ml 7©^ 



iDonati et aLld 19971) measured a vsini of 12 + 2 km s which 
suggests either that HD 10 4237 is viewed clo se to pole-on or 
that it is a moderate rotator iGradv et al] (l2004l) found indeed an 
inclination of the stellar rotation a xis of 18°;^}p which implies 
an equatorial velocity of 38 km s~' . iBohm et al.l (1200 6^ detected 
a modulation of the Wa line profile with a period of 100 ± 5 hrs 
and suggested a rotational origin, yielding a similar incUnation 
of 23°!g. This value of P.^t was confirmed bv lTesta et aP d2008l) 
in a X-ray study of the HD 104237 system. 



Concerning magnetism. lDonati et al.l d 19971) detected a weak 
field of about 50 G in HD 104237, but a later survey by 
IWade etalJd200l could not confirm it at this stage. 

As can be seen from all the different contributions cited here- 
above, no in-depth spectral type determination has been done as 
of today: indeed, all determinations were based either on pho- 
tometry o r on low-resolution spectroscopy (except the rough es- 
timate bv lDonati et alJ[T997[) . and a thorough redetermination of 
fundamental parameters based on high-resolution spectroscopy 
is required for an independent conclusion. This is one of the ma- 
jor goals of the work presented in this article in Sect. [6[ and 
a basic requireme nt for a subsequent asteroseismic modeling 
dFumel et al.ll201 ih . but also in order to precise the PMS insta- 
bility strip features and boundaries. Knowing these parameters 
will eventually allow us to better understand the origin of the 
activity observed in the Herbig Ae/Be stars and test the current 
theoretical evolutionnary models for PMS stars. 
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Table 1. Stellar parameters of HD 104237. 2'"' column: P corresponds to the stellar parameters of the primary and S to those of 
the secondary. Values in italics (log ^ and R/Rq): computed from universal gravitation and Stefan-Boltzmann laws. References: 
[11 [Hu et al. (1989), [2] Hu et al. (1991), [3] Brown etal. (1997), [4] van den Ancker et al. (1997), [5] Donati et al. (1997), Jj6] 
Ivan den Ancker et al. (1998), [7] Gradv et al. (2004), [81 iBohm et al.1 (.2004). [91 ILuhmani (|2004 . [101 lAcke & WaelkensI (l2004h . 
[Ill Lvo et al. (2008), [121 Bohm etal. (2006). M ass and age s are d ete rmined using the mod el i sochrones and isomass of: 
llben & Renzini (1984), C'' Palla & Stahler ( 199i, ''^' ISiess et al.1 (l2000l) . ^ iPalla & Stahie3(IT999l) . '''' iPalla & Stahie3(l200l . 
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2.2. Previous pulsational analysis 

The primary component of the HD 104237 system has revealed 
being a pulsating Herbig Ae star From spectropol arimetric ob- 
servations of active stars including HD 104237, iDonati et al.l 
(1199 7) detected for the first time radial velocity variations in 
equivalent photospheric LSD (cf. Sect.[3]l Stokes 1 spectra of this 
star and determined a period of 37.3 min. This short period was 
incompatible with a rotational modulation or a secondary com- 
panion and was attributed to stellar pulsations. This pulsational 
variability, typical of 6 Scuti-ty pe vari ables, was confir med with 
photometric measurements bv iKurtz & Miiller (1 19991) who de- 
tected two frequencies. They computed the ratio of radial veloc- 
ity to V light amplitude and found values typical for a S Scuti 
star. Moreover, they calculated a value for the pulsation constant 
Q lower than the typical values given by the fundamental, the 
first and the second overtone in 6 Scuti stars, indicating pulsa- 
tions of higher overtone. 

iBohm et al.1 (|2004 obtained in 1999 and 2000 an extensive 
high-resolution (R«s35000) spectroscopic data set of HD 104237 
at the 1.9m Radcliffe telescope of the SAAO (see Sect.[3ll. The 
resulting high quality radial velocity curve allowed them to de- 
tect for the first time by spectroscopic means multi-periodic os- 
cillations in a PMS star (8 frequencies in the 1999 data set and 
5 in the 2000 data set). Five frenquencies of 1999 were identi- 
fied with frequencies of 2000, between 28.50 and 35.60 d"'. No 
pulsation was detected in HD 104237b in the night around the 
periastron, when both components are well separated in veloc- 
ity, a result which is not astonishing given the position of the star 
in the HR diagram. The frequency differences of the 8, respec- 
tively 5 frequencies detected by Boh m et al.l (12004 ) were smaller 
thant the estimated large separation (Avq 40 - 45 yuHz), lead- 
ing the authors to conclude on the nonradial nature of at least 
some of the pulsation modes detected in their work. More gen- 
erally, no regular frequency pattern has been clearly identified 
in the pulsational behaviour of HD 104237. iBShm et all (t2004 
pointed out (based on the frequency separation) the potential ex- 
istence of 1 to 3 radial pulsation modes among the 5 frequencies 
observed both in 1999 and 2000, but at this time no clear mode 
identification in HD 104237 could have been performed for none 
of the modes. Since this first article was restricted to the analysis 
of radial velocity variations, we wanted to go a step further, and 



reanalysed the spectral line profiles themselves in the same data 
set; the result of this study is presented thereafter. 

The previously determined pulsation characteristics are sum- 
marized in Table 2\ The three last frequencies of 1999 from 
iBohm et al.l ^mife to /s) are included, although these au- 
thors considered them as les s certain . 

Using the results of Bohm et al.1 (|2004 . iDupret et all (l2006l 
I2OO7I) found in a preliminary seismic modeling of HD 104237 
that it is difficult to properly fit the observed frequencies 
with theoretical frequencies computed from mode l s with Teff 
and log(L/L o) presented in Ivan den Ancker et al.l (11998 ) and 
iGradv et al] (12004,) . An important result was their conclusion 
that HD 1043237 is in fact not a typical 6 Scuti type pulsator 
Indeed, the p modes order corresponding to the observed fre- 
quencies are too high to be excited by the standard excitation 
mechanism of 6 Scuti type variables, namely a A--driving in the 
He 11 partial ionization zone. Such a modeling requires a pre- 
cise knowledge of the stellar fundamental parameters which en- 
ter the evolution and oscillation codes and a good idea of { and 
m for each observed mode. However, both are poorly known 
or not known at all in the case of HD 104237. In order to go 
further in this seismic study and to better constrain the inter- 
nal structure of HD 104237, we will propose a more thorough 
and complete modeling of HD 104237 in a forthcoming article 
(Fumel et al. 2011), using the accurate stellar parameters deter- 
mination as well as the pulsational information described in the 
present article. 



3. SAAO spectroscopic observations and data 
reduction 

The analysis presented in this artic l e is b ased on parts of the 
data set described in iBohm etal] (|2004 . namely the high- 
resolution spectroscopic observations obtained in April 1999 
and April 2000 at the 1 .9 m SAAO (South African Astronomical 
Observatory) Radcliff telescope; all data were obtained with the 
GIRAFFE fiber-fed echelle spectrograph with a resolving power 
of about/; ^ 35000. 

HD 104237 was monitored in quasi-continuous mode 
throughout all the nights of the 1999 and 2000 runs (resp. 7 
and 14 nights) in order to obtain high-resolution spectroscopic 
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Table 2. Previous pul sational analysis res ults for HD 104237. 
Col. 5 : references: (1) 'Ponati et al.l (Il997h . (2) iKurtz & Miilleil 
(Hm, (3) Bohm et al. (2004). 



Year of 


Frequency 


Period 


Amplittide 


Kei. 


observation 


(d ) 


(min) 






1993 


/i = 37.3 


38.6 ± 1 


0.65 km S-' 


(1) 


1995 


/, = 39.6 


36.4 ± 1 


0.65 km s"' 


(1) 


1998 


/i = 33.29 


43.3 


11.2 ± 0.5 mmag 


(2) 




/2 = 36.61 


39.3 


3.4 ± 0.5 mmag 




1999 


/i = 33.289 


43.257 


1.320 km S-' 


(3) 




/z = 35.606 


40.443 


0.474 km S-' 






/3 = 28.503 


50.521 


0.195 km s"' 






/4 = 30.954 


46.521 


0.139 kms"' 






/s = 33.862 


42.525 


0.099 km 






/s = 32.616 


44.150 


0.105 km S-' 






'fi = 34.88 


41.28 


0.1 kms"' 






fs = 35.28 


40.82 


0.05 km S-' 




2000 


fi = 35.609 


40.439 


0.328 km S-' 


(3) 




/z = 33.283 


43.265 


0.258 km S-' 






/3 = 31.012 


46.434 


0.177 km S-' 






/4 = 28.521 


50.489 


0.165 km S-' 






A = 32.375 


44.479 


0.113 kms-' 





time-series necessary for an asteroseismic analysis. The covered 
wavelength domain usually was 426-688 nm (spread over 50 or- 
ders) in 1999 and 436-688 nm (spread over 47 orders) in 2000. 
In order to sample properly the previously known main pulsa- 
tion period of about 40 minutes, subexposure times of 5 minutes 
were needed, yielding about 8 spectra per pulsation period, the 
brightness of HD 104237 enabling such short exposure times. 
A standard caUbration strategy was applied. Values of Signal to 
Noise Ratios (SNR) per pixel at 550 nm ranged from 30 to 120 in 
1999 and from 30 to 1 10 in 2000, with typical values of 50-70. 
Most of the data reduction was carried out following standard 
reduction pro cedures using the " ESPRIT" spectroscopic reduc- 
tion package (iDonati et al.lll997h . Moreover, as the analysis of 
pulsations by high-resolution spectroscopy requires significant 
S NR, the least-square deconvolution (LSD) method as described 
in iDonati et al.l d 19971) was applied. This method makes use of 
the multiplex gain of more than 500 lines present in the spec- 
trum of this star, providing a high SNR equivalent photospheric 
profile. More precisely, it assumes that the local profile of all 
selected spectral lines is similar in shape (with different individ- 
ual central depths), the line parameters being extracted from an 
appropriate catalogue, corresponding as close as possible to the 
fundamental parameters of the star of interest known at this stage 
of d ata reduction as A4-type (based on the results announced in 
IDonati et al. 1997). In a subsequent step, the more than 100 nar- 
row telluric vapor lines present in the spectra yielded the nec- 
essary multiplex information for a high precision radial velocity 
correction. All spectra were shifted to the heliocentric rest wave- 
length, and calibrated to a final intrinsic precision estimated to be 
around 100 ms ' for this data set. In addition, HJD (Heliocentric 
JuUan Date) have been calculated for each LSD profile. Details 
of observations and data reduction are thoroughly described in 
[Bohm et al. (2004). 

The detection and identification of non-radial pulsations de- 
scribed in Sect. |4]has been carried out on the data collected dur- 
ing the SAAO observing run in 1999, where we especially con- 
centrated our study on the data of the nights from 22"'' to 26* 
of April, these nights having yielded the longest individual time 
coverages. During the night of April 25* 156 high-resolution 
spectra were obtained during 10.9 hrs. 



The fundamental parameter determination described in Sect. 
|6]has been performed on data collected during the night of 12* 
of April 2000 at SAAO; on this d ate, the eccentric binary orbit 
with a period of 19.859 days (see iBohm et al.l 12004 ) separates 
during periastron both spectra with a relative shift of 60 kms 
thus avoiding at maximum pollution of the primary spectrum by 
the much fainter secondary companion HD 104237b. 

The log of the observations we worked with are summarized 
in Tabled 



Table 3. Log of the observations we worked with. 



Date 


Time series 


Number of 


SNR (pixel-') 




duration (hrs) 


spectra 


at 550 nm 


22 Apr. 1999 


8.2 


74 


55-120 


23 Apr. 1999 


10.2 


91 


40-90 


24 Apr. 1999 


10.8 


136 


30-100 


25 Apr. 1999 


10.9 


156 


30-60 


26 Apr. 1999 


7.7 


58 


30-80 


12 Apr. 2000 


3.6 


35 


40-60 



4. Detection and identification of non-radial 
pulsations 

Stellar pulsations are characterized by their frequencies and cor- 
responding modes, specified by three quantum numbers: the or- 
der n, related to the number of radial nodes, the degree of the 
mode £, indicating the total number of surface nodes, and the az- 
imuthal order of the mode m such that m - -£, -£+ 1, ...,£- 
where \m\ coiTesponds to the number of meridian nodes, i.e. 
showing nodal lines perpendicular to the stellar equator. Radial 
modes are characterized hy £ = m = and non-radial modes 
by ^ 0. Because of cancellation eff'ects, photometry or ra- 
dial velocity measurements are not able to access higher degrees 
of non-radial pulsations with typically £ > 3 (iKennellv et al.l 
11996). The only direct way to investigate modes of higher de- 
grees consists in analysing line-profile variations of rotationally- 
broadened photospheric profiles, where the 3D velocity field of 
the oscillating stellar surface induce s profile variations due to 
Doppler effect (IVogt & Penrodlll983h . The maximum attainable 
degree £ is increasing with the ratio v sin / to resolved spectral 
element (Av - cjR, with R being the resolution of the spectro- 
graph). This also implies that fast rotators seen almost pole-on 
only allow mode detection limited to low degrees, due to a lack 
of line profile width. 

As explained in previous Section, we concentrated our study 
on the data of the 22"''-26* of April 1999. The initial work by 
iBohm et alJl2004l did not reveal any non-radial pulsation at first 
sight in the line profiles and these authors therefore decided to 
concentrate on radial velocity variations, also due to the fact that 
the star is a moderate rotator, seen 3/4 pole on and exhibiting a 
small y sin/ of only 12+2 kms However, line profiles of this 
star are broadened in the wings by some additional agent, lead- 
ing to a total width of almost 50kms-', the double of the pure 
rotational broadening. The resolution of the GIRAFFE spectro- 
graph at SAAO being of the order of R 35000, a resolved ele- 
ment corresponds roughly to 8 km s"' . This implies that a search 
for non radial pulsations (NRP) in this instrumental context is 
limited in this star to degrees £ of the order of 2 to 3, higher 
degrees being not observable due to a lack of rotational veloc- 
ity resolution. In the present work we decided nevertheless to go 
one step further, and to reanalyse in depth the profile variations 
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themselves, not limiting the study to radial velocity variations. 
To do so, and in order to search for faint traces of non-radial 
pulsations in the primary component, each LSD profile has been 
corrected for its radial velocity. This correction comprises orbital 
motion, but also centroid shifts mainly due to radial pulsations. 
After recentering, a mean LSD profile, averaged over the whole 
run, has been subtracted from every profile in order to search 
for very faint signatures. This allowed us to detect line profile 
variations clearly indicative of non-radial stellar pulsations, in- 
ducing differential amplitudes of the order of 1 .5% of the contin- 
uum. Also, it is obvious from the differential dynamical spectra 
of FiglT]that beating between nearby frequencies occur, which 
leads to the conclusion that probably more than one NRP must 
be acting (the night of April 26* is not represented due to shorter 
coverage). The sampling of the LSD profiles in dispersion di- 
rection satisfies the Nyquist criteria, i.e. one resolved element 
projects on 2 velocity bins of 4kms"'each. Fig|2]shows the dif- 
ferential amplitude variations at line center, extracted by taking 
the median value of the 3 central velocity bins of each spectrum. 

HDl 04237 SAAO 1999 observations 




AHJD 



Fig. 2. Amplitude variation of the central velocity bins of the 
residual LSD profile of the nights of Apr 22"'' to 26* 1999. 
Superimposed is a fit of the data corresponding to the result of 
the frequency analysis as presented in Tab. |4] 



the frequencies, amplitudes and phases of multi-periodic signals 
in not-equally spaced data sets. Both process fit simultaneously 
all the frequencies, amplitudes and phases detected so far when 
a new peak is found in the DFT spectrum of the time-serie. 
Unlike Period04 which provides a flexible interface to perform 
multiple-frequency fits, SigSpec is fully automatic. In addition, 
it includes a rigorous statistical treatment of how to compute the 
significance level of these peaks (with respect to white noise) 
and provides directly a significance parameters, the spectral sig- 
nificance sig. In order to avoid at maximum low frequency night- 
to-night variations, we corrected our central velocity bin time 
series by subtracting a nightly average. Still, by doing so, some 
low frequencies persisted as can be seen later on. 

The results of the Period04 analysis are presented in Table |4] 
based on th e amplitude varia tions as shown in Fig. |2] As 
described in iBohm et al.ll2004l the uncertainty in the extracted 
frequency value can be estimated in differe nt ways. The mos t 
conservative approach has been proposed by dRipepi et al.ll2003h 
who suggest to estimate the error by measuring the FWHM 
of the main lobe of the spectral window function, correspond- 
ing here to 0.31 d"', while ( Kurtz & Miiller 1999) use as an 
estimator the time basis ISJ of the observing run, yielding 
for the SAAO 19 99 run a much sma ller error estimate of 1/(4 
Ar) = 0.04d iBreger et afl (Il993h studied empirically the 
possibility of a peak in the periodogram being a true signal of 
pulsations with respect t o the n oise level; this work was refined 
later by Kuschnig et al.l (Il997h : the significance of a peak in 
the amplitude periodogram exceeding 4.0 times the mean noise 
amplitude level after prewhitening of all local frequencies has 
a 99.9% probability to be due to stellar pulsations (99.0% for a 
ratio of 3.6, 90.0% for a ratio of 3.2). Fig.[3]shows the amplitude 
power spectrum including all strong peaks around 29 - 36 d 
the noise levels having been determined after prewithening of 
the 9 first frequencies, while Fig. |4] shows the fact that even F9 
has a signal to noise level exceeding 4.0. We arbitrarily decided 
to stop the iterative procedure for amplitude values below 1%. 
The verification of the results with SigSpec revealed exactly the 
same frequencies and amplitudes, all detected frequencies Fi 
to Fio having a significance level between 66.5 and 6.3, respec- 
tively. No combined frequency is present in the power spectrum, 
following the analysis with the Combine program by Ree genI 
(l2007h . We don't attribute much importance to the frequencies 
below approximately 5d"', since differential calibration issues 
from one night to the other might persist, i.e. this excludes F4, 
F7, Fs and Fm- By applying the more conservative error estimate 
by (iRipepi et al.l 12003;) . all remaining frequencies ("F") could 
be identified with the frequencies ("f ) determined by ra dial 
velocity studies on the same data set in iBohm et al.1 (I200 4'). F3 
being slightly outside the error bar, and the 1.0d~'difference 
between F9 and f3 indicates that we are most likely in presence 
of an alias of the same frequency. It can be noticed, that the 
order of the first two dominant frequencies is inverted depending 
if radial velocity or central depth variations are analysed. 



4.1. Central profile analysis with SigSpec and Perlod04 

We choose to carry out a frequency analysis of o ur LSD residual 
profil es time-series using the Period04 package (Lenz & Brege^ 
l2005b in parallel with the Sigspec package (Reegen 2007). They 
both use, in iterative process, combination of Discrete Fourier 
Transform (DFT) and least-squares fitting algorithms to extract 



4.2. Mode Identification with the Fourier 2D method 

A direct method of analysing the non-radial pulsation modes 
present in the LSD-spectra time series consists in applying the 
Fouri e r 2D method o n line-profile variation s dKennellv et al.l 
1 19931: iKennellvl Il994t iKennellv et aP Il996l) . This technique 
analyses the complex pattern present in the line profiles by 
computing a 2-dimensional Fourier transform in both time and 
Doppler space. To do this, an interpolation of each profile on a 
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Fig. 1. Line profile variations due to the non-radial pulsations of HD 104237. The deviations from the mean intensity are displayed 
for four different nights (top left: April 22"'' 1999, bottom left: 23"^ 1999, top right: April 24*, bottom right: 25* 1999). 



grid representing stellar longitudes is performed, transforming 
velocities across the line profile into longitudes on the stellar 
equator using the relation Av = v sin / sin (p. Av is here the veloc- 
ity position within the LSD profile with respect to the rest wave- 
length of the star and (/> is the stellar longitude angle of the star in 



spherical coordinates. In the resulting two-dimensional Fourier 
spectrum, the temporal frequencies are related here mainly to 
the frequencies of oscillation, while the apparent azimuthal or- 
der m is related to the structure of the modes present at the stellar 
surface, without being identical to the usual azimuthal order m. 



7 



A. Fumel and T. Bohm: Spectroscopic monitoring of the Herbig Ae star HD 104237 



Table 4. Frequencies determined with Period04. The different 
columns are: (1) number, (2) frequency in d^'and (3) in /zHz. 
(4) Amplitude in (F/Fc) (flux variation with respect to the aver- 
age pr ofile), (5) identification with frequencies as of lBohm et alJ 
(12004 . and absolute shift in d 



# 


FpQ4 


PpQ4 


ApQ4 


fB04 (l<5f 1) 




(d-') 


0"Hz) 


(F/F,) 




F, 


35.60 


412.04 


0.0089 


fj (0.01) 


Fo 


33.74 


390.51 


0.0033 


f, (0.12) 


F3 


32.25 


373.26 


0.0030 


f6 (0.36) 


F4 


4.47 


51.74 


0.0024 




F5 


31.08 


359.72 


0.0028 


f4(0.12) 


F6 


34.00 


393.52 


0.0019 


f5(0.14) 


F7 


1.94 


22.45 


0.0018 




Fg 


2.72 


31.48 


0.0014 




F, 


29.48 


341.20 


0.0011 


f3 (1.0) 


Fio 


5.17 


59.84 


0.0010 





HD 104237 - F1 SAAO 1999 




Frequency (c/d) 



Fig. 3. Periodogram of the line center variation as shown in 
Fig|2] without prewithening. The large buUc of frequencies 
around 29 - 36 d^'can be seen, as well as the 3.2, 3.6 and 4.0 
mean-ampitude level of the noise determined by prewithening 
Fi to F9. 

HD 104237 - F9 SAAO 1999 




Frequency (c/d) 



Fig. 4. Same figure as Fig. [3] but prewithened by Fi to Fg. Only 
Fg is present, well above the 99.9% confidence level. 



The original work bv lKennellvl(ll994 showed in fact that appar- 
ent \m\ scales as i +2 for values close to zero, as ^ -i-l for val- 
ues lower than 10 and as £ for values above 10. We performed 
this F2D computation on the night of April 22"'' to 25* 1999. 
A weighted c ombination of the n ightly F2D spectra can be seen 
in Fig. l5] (see iBohm et al.ll2009l for details of the method) and 



our results confirm the presence of low-degree non-radial pulsa- 
tions: a peak at the dominant frequency Fi is clearly seen with 
an apparent m comprised between 2 and 4 indicating a degree f 
comp rised between an d 2. Moreover, as can be seen e.g. in Tab. 
2.3 of lKennellvl (fT994l) (or in simulation with available online- 
NRP simulators), the symmetric pattern of the dynamical resid- 
ual profile as seen in Fig[T]indicates for Fi an azimuthal order m 
= ± 1, implying a most likely C value of 1 or 2. 

The frequency peak in Fig. |5]between 3 and 5d 'is again 
not taken into account, since nightly calibration shifts can easily 
introduce such variations. 



2D Fourier analysis HD1 04237 SAAO 1999 




Fig. 5. Combined F2D analysis of HD 104237 for the nights of 
Apr. 22"'' to 25* 1999. 



A complementa ry analysis on the whole da ta set using 
the FAMIAS code (Zi ma et all l2006t IZimal l2008l) is planned, 
but the highly asymmetric mean photospheric profile seen in 
HD 104237 adds complexity in the data analysis and cannot in 
the actual state of the code be taken into account. 



5. Optimized spectral continuum determination 

In addition to the best available mode identification, a future 
modeling of the stellar oscillations of HD 104237 requires the 
precise redetermination of its fundamental stellar parameters. 
Since we decided to carry out this analysis from reliably mea- 
sured equivalent width of selected absorption lines (see Sect.|6]l, 
a perfect continuum normalization revealed to be essential, par- 
ticularly in order to significantly increase the SNR by summing 
all the spectra of the night. An automatic code was therefore de- 
veloped, the case of HD 104237 being rather complicated due 
to its spectroscopic binarity: the pollution of the primary (P) 
spectrum by the secondary (S), but also the presence of emis- 
sion components and variability observed in numerous absorp- 
tion lines require particular attention in the normalization proce- 
dure. The data set obtained at SAAO in 2000 (see Boh m et al.l 
2004) contained one night of observations very close to the pe- 
riastron of the binary (12* of April). We decided therefore to 
concentrate on the data of this particular night in which both 
spectra are well separated in velocity, making easier to distin- 
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guish P from S. As we worked on echelle spectra, we carried out 
an individual normalization of each order, performed in 4 mainly 
automated main steps. 

As a first step, we created a comb including efi'ective con- 
tinuum locations and eliminating unsuited others areas. To do 
so, we extracted 2 lists of absorption lines (Teff (P) = 8500 K 
and reff(S) - 4750 K, \ogg - 4.0, Vmicro - 2 km s"' and solar 
abundan ces) from the Vienna Atomic Line Data b ase (VALD; 
Piskuno vet all 119951: iRvabchikova et all 119971: iKupka et all 
1999i.i2000l) on which we apphed: 



i) a selective criterion on the maximum depth accepted for spec- 
tral lines to be considered as continuum (< 1% of the contin- 
uum), taking into account the relative luminosities of both com- 
ponents as a function of wavelengths, assuming two blackbodies 
with temperatures equal to 8500 K and 4750 K, 

ii) a criterion on the atypically large profile width of +2.5 v sin / 
excluding areas with spectral lines, 

iii) the elimination of specific areas whith emission, variable, tel- 
luric and the broad Balmer lines (the latter extending over sev- 
eral spectral orders and therefore not reliably to be normalized), 

iv) a shift of P and S wavelength grids by respective orbital ve- 
locities at the given observation time (all orbital p arameters of 
the sp ectroscopic binary have been determined in iBohm etal] 
12004 . 

After multiplication of the resulting comb with the observed 
spectrum, only wavelengths ranges corresponding to the remain- 
ing continuum areas were kept, all other areas were rejected. The 
next step consisted in fitting a polynomial of degree 4 through 
the remaining points in each order and a division of each order 
of the observed spectrum by this optimized polynomial. Finally, 
in a last step the software procedes to a weighted concatenation 
of the normalized orders, the overlap between orders i and i-i-1 
were taken into account through the following equation : 



Foveilap(^) 



Fi(A) X SNRf(i) + Fm(A) X SNRf^j(^) 
SNRf(^) + SNR2^j(i) 



(1) 



where Foveriap(^) is the total flux resulting from the addition of 
two overlapped parts of consecutive orders with normalized flux 
F{A), weighted by their corresponding local SNR(/l), resulting 
from a polynomial fit of the respective SNR curve within the 
overlap, the local SNR being part of the reduced spectra files. 

Since Balmer line profiles span several orders, they were not 
taken into account in this normalization process, i.e. area from 
4816 to 4929 A and from 6475 to 6645 A remained unnormal- 
ized. During periastron both binary components show significant 
radial velocity shifts within one night. In order to sum up the 
large number of individual spectra during that night, we decided 
to center each individual spectrum to the rest wavelength of the 
primary component. Summation was thereafter done in the fol- 
lowing way: 



F,ot(^) = 



i:]Li Fj(^) X SNR2 



iSNR2 



(2) 



where Ftot(/l) is, at each wavelength, the summed normalized 
spectrum, n is the number of spectra of the night, F{A) the in- 
dividual flux, SNR is the mean SNR of each spectrum. After 
summation, the gain in SNR of our spectrum was of ^ 6 since 
the SNR per resolved element at 550 nm increased from 64 to 
372 (from 45 to 263 per pixel; see Fig.|6ll. 
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Fig. 6. Individual (thin line) and summed (thick line) normal- 
ized spectrum of the periastron night around one of the selected 
lines for our study (Fel 5415. 1920 A). Synthetic spectra with 
Teff =8500 K (dotted fine) and T^g =4750 K (dashed line) have 
been overplotted (log ^=4.0, v^i^o = 2kms"'solar abundances 
for both spectra), showing the location of the individual lines. 
The secondary spectrum has been corrected by the local lumi- 
nosity ratio. The figure shows that the line selection process de- 
scribed in Sect. |6] works well. However, it can be noticed that 
either fundamental parameters and/or local luminosity ratios are 
not yet perfectly tuned at this stage. 



6. Fundamental parameters determination 

An accurate and reliable knowledge of stellar fundamental 
parameters, mainly effective temperature T^ff, surface grav- 
ity \ogg and chemical abundances log A (defined as follows: 
logAE = log(^) for an element E) is crucial to constrain the 
stellar atmosphere model, which is the basic ingredient of a 
forthcoming asteroseismic modeling. As shown in Sect. |2l for 
HD 104237 large discrepancies persist between parameters an- 
nounced in the literature. A precise redetermination of its funda- 
mental parameters is therefore crucial. 



6.1. Method 

Since Teff, log^, chemical abundances (logA) and microturbu- 
lent velocity Vmicio are interdependant and affect simultaneously 
the depth, width and shape of spectral absorption lines, they have 
therefore to be determined simultaneously. Numerous methods 
exist for such an analysis, but very few of them are suitable in 
the particular case of HD 104237. Due to the non standard shape 
of HAeBe spectral energy distributions (IR excesses, UV deple- 
tion), classical photometric determination of T^s have to be con- 
sidered with caution. Another frequently used method, the stellar 
parameter determination using Balmer lines have to be excluded 
because of core emission in the line and/or P Cygni profiles, but 
also on a technical side due to normalization issues in case of 
high dispersion echelle spec tra. Techniques ba sed on lines ratio 
for determining Teff (see e.g. lSousa etal ] |20T(I require non pol- 
luted lines and were rejected because of strict selection criteria 
(see Subsect. I6.1.3I I. which excluded all workable lines in our 
case. 
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Table 5. Characteristics of absorption lines selected for fundamental parameters determination. Column are: (1) central wavelengh, 
(2) ion, (3) excitation potential (eV), (4) log of statistical weight x oscillator strength of the line transition, (5) local mean Signal to 
Noise, (6) observed equivalent width, (7) absolute and (8) relative uncertainty on the equivalent width. 





el.+ion. 




log^f/ 








rr f'PW A 
(TjelttlWobs; K'O) 


4485.675 


Fel 


3.686 


-1.020 


111 


16.1 


2.6 


15.9 


4602.941 


Pel 


1.485 


-2.209 


245 


39.3 


4.0 


10.1 


5049.819 


Fel 


2.279 


-1.355 


268 


51.9 


3.9 


7.6 


5132.669 


Fell 


2.807 


-4.094 


275 


28.0 


4.0 


14.5 


5393.167 


Fel 


3.241 


-0.715 


289 


46.7 


3.1 


6.6 


5415.192 


Fel 


4.386 


0.642 


285 


97.9 


4.4 


4.5 


5434.523 


Fel 


1.011 


-2.122 


273 


51.9 


3.3 


6.3 


5576.089 


Fel 


3.430 


-1.000 


272 


35.2 


3.5 


9.8 


6084.111 


Fell 


3.199 


-3.881 


245 


18.2 


4.0 


21.8 


6400.000 


Fel 


3.602 


-0.290 


238 


51.3 


4.7 


9.2 



After consideration, we favoured an equivalent width (EW) 
analysis with respect to a direct spectrum fitting method, be- 
cause of the existence of additional broadening agents and strong 
asymmetries in the line profiles, whose origin is unknown and 
which are therefore difficult to reproduce, and of numerous ac- 
tivity features in the whole spectrum that complicate direct line 
profile fitting. Our method consist in a comparison of EW of a 
set of selected lines, between the observed combined spectrum 
of the periastron night of HD 104237 and a fine grid of synthetic 
spectra in the 3 fundamental parameter dimensions, namely T^s, 
log g and chemical abundance of a given chemical element (re- 
quired to work out reliable values of Teff and log^). With more 
than 2000 spectral lines featured in a 8500 K synthetic spectra, 
Fe turned out to be the most suitable chemical element to estab- 
lish our EW analysis. The high activity level of this HAe star 
(emission, variability, line asymmetry) and its binarity make a 
reliable determination of fundamental parameters in this case 
very challenging. 



6.1 .1 . Correction of the secondary spectral contribution 

Since our main goal was to determine the fundamental param- 
eters of the primary component, we had to correct for the pol- 
luting spectroscopic contribution of the secondary component, 
HD 104237b. Indeed, although the secondary is significantly 
fainter than the primary component, it is also much cooler and 
does exhibit strong metallic absorption lines, leading to system- 
atic contamination of primary lines and even lowering in some 
parts of the spectrum the continuum level. To free the observed 
spectrum from this effect, a synthetic spectrum corresponding 
to the secondary component was substracted from the observed 
spectrum. In order to minimise additional errors from an incor- 
re ct model of seconda ry, we used stellar parameters determined 
bv lBohmet all ( 12004 based on cinematics of the binary move- 
ment, adopting solar chemical abundances, a log^ value of 4.0 
and a Vmicro of 2 kms"': several synthetic spectra, with Tgs{S>) 
comprised between 4250 K and 5250 K (by step of 250 K) were 
therefore tested. Since the observed normalized spectra had been 
centered on the rest wavelength of the primary component, sec- 
ondary synthetic spectra were shifted by the respective orbital 
velocity displacement. We called a{A) a monochromatic lumi- 
nosity ratio with a(A) - B^{X)I Bs,{A), where B^(A) and B^{A) are 
respectively the Planck laws of the primary and the secondary. 
The most likely secondary model was determined by computing 
and minimizing a quantity between the observed spectrum 
and the expected secondary contribution, yielding a T^s value of 
4500 K for the secondary. A spectrum corresponding to this tem- 



perature (Fs(/l) flux) was therefore subtracted from the observed 
spectra (Fobs (/I) flux), with a contribution based on the equation: 

Fp(i) + Fs(^)/a(^) 



Fob,s(^) = 



(3) 



1 + 1 la(A) 

in order to obtain a primary spectrum (Fp(/l) flux) cleaned up 
from the secondary contribution. Obviously, systematic uncer- 
tainties of the monochromatic luminosity ratio and the intrinsic 
stellar parameter selection of the secondary do persist and are 
difficult to fully be taken into account. 

6.1 .2. Construction of a 3D-grid of syntlietic spectra 

In order to precisely redetermine the fundamental parameters of 
our star we needed to compare the observed high SNR spectrum 
to synthetic spectra. Therefore, we constructed a 3-dimensional 
grid (Teff, \ogg, logApe) of synthetic spectra based on line cat- 
alogs provided by the VALD database (see Se ct|5l) as input of 
the SYNTH3 spectrum synthesis code dKochuk hov 2007), each 
spectrum covering a wavelength range from 4400 A to 7000 A. 
The S YNTH3 code assumes Local Thermodynamic Equilibrium 
(LTE) and a plane-parallel hydrostatic stellar model atmosphere. 
For early type stars input models are provided by the well known 
Kurucz program suite. We selected solar abundance atmosphere 
models from the grid, after having checked that the effect of an 
input model containing slightly different than solar metallicities 
is of second order on the calculated synthetic spectrum. 

As a starting point, we adopted parameters derived 
by previous authors to define the boundaries of the grid: 
7500 K < Teff < 9500 K (by steps of 250 K) and 3.5 < log g < 4.5 
(by steps of 0.5). Since it is not uncommon to find HAeBe stars 
exhibiting in their spectra simultaneously the presence of some 
elements in suprasolar abundance, whereas other ones are sub- 
solar (see e.g. Acke & Waelkens 2004), without particular trend 
in terms of global metallicity, we decided to modify in this 
study only the Fe abundance in our VALD requests. It should 
be noticed at this stage, that studying the iron abundancy was 
a necessary by-product in order to constrain effective tempera- 
ture, as we will see subsequently, log Ape values comprised be- 
tween -4.84 and -4.00 seemed to be reasonable, i.e. [Fe/H] 
values comprised between -0.3 and -1-0.54 with respect to the 
solar Fe abundance given in VALD, which adopted logApeg = 
-4.54. As a reminder, [Fe/H]=log(^)^ - log(^)^. Typical 

values of Vmido found in the literature were of 2-3 km s~'(see 
e.g. Acke & Waelkens 2004; Guimaraes et al. 2006; Cata la et al.1 
2007]), we therefore fixed a value of 2 km s 'for the microturbu- 
lent velocity HD 104237. A wavelength by wavelength quadratic 



10 



A. Fumel and T. Bohm: Spectroscopic monitoring of the Herbig Ae star HD 104237 



interpolation was subsequently done between synthetic spectra 
in order to obtain a finer spectral grid with a 25 K step in Teff, a 
0.1 step in log^ and a 0.01 step in logApe. 



6.1.3. Photospheric line selection 

Photospheric iron absorption lines used in our study were se- 
lected using the following criteria: 



- exist for every (T^fi, \ogg, log Ape) triplet in order to be in- 
tercomparable; 

- exclude blends with lines from other chemical elements 
(than Fe), i.e. within a range of +2.3 v sin i/c from the central 
wavelength (the extraction width having been adoped after 
visually inspecting the spectra and taking into account the 
broad wings of the photospheric lines of our star); 

- avoid pollution by lines from the secondary spectrum, even 
if a best-fit secondary spectrum had been removed, in order 
to minimise the impact of the potential error on reff(S). To do 
so, secondary lines with a relative contribution deeper than 
10% of the central depth led to a rejection of the primary 
spectrum line. 

- show no or little variability in time (especially over short 
periods of time); 

- do not exhibit at any time an emission or other than normal 
(for HD 104237) asymmetry; 

- be strong enough with respect to local noise, in order to min- 
imize relative uncertainties on EW and to be able to neglect 
residual contribution from very faint lines of primary and 
secondary spectra (VALD requests require a depth thresh- 
old and provide catalog of lines deeper than this thresh- 
old value), but not to strong to avoid saturation effects. 
As a result, we kept weak and moderratly strong lines 
(15mA<EWobs < 100 mA); 

- show exploitable nearby continuum, e.g. with no emission 
component in their immediate vinicity (after automatic se- 
lection, a individual visual check line by line was done); 

- cover a variety of excitation potentials, ionization stages and 
wavelengths; 

- having ascertained atomic pa rameters, in partic ular 
log;?/ values (comparison with iFuhr & Wiesd l2006l and 
lMel6ndez & Barbuv 2009) . 



This rigorous selection resulted in a set of only 10 spectral lines 
(8 Fe 1 lines and 2 Fe 11 lines) which are listed in Table |5] 



6.1.4. EW measurements 

EW measurements of our selected lines were carried out using a 
trapezoidal integration. For consistency, the identical extraction 
limits were used for the integration of a given line from both the 
observed spectrum and the synthetic spectra. 

To extract errors on EW values, we assumed that they 
mainly came from local averaged photon noise and we ne- 
glected uncertainties in integration limits or arising from 
residual continuum normalization errors. Multiplying the lo- 
cal SNR value by the line integration width Av of each line, 
we then obtained an estimation of absolute uncertainty in 
EW, crabs(EWobs), and its corresponding relative uncertainty 
cri-ei(EWobs) = crabs(EWobs)/EWobs (see Table|5]l. 



6.2. Teff, logg and Fe abundance determination 

In order to determine T^g, logg and log Ape of HD 104237, it 
was necessary to compute the degree of similarity between our 
observed spectrum and all the synthetic spectra of our 3D-grid 
by means of the 10 previously measured EW; a direct synthetic 
spectrum fitting technique could not be applied due to highly 
atypical photospheric line profiles (asymmetry, broad wings). 

A maximum likelihood estimation was therefore done by 
minimizing the following reduced merit function 5 red: 



J=l ^ 



/ E Wobs, j - EWsynth, j ( T^eff , log g , log Ape) \^ 



(4) 



where N-k is the number of degrees of freedom, N being the 
number of selected lines and k being the number of param- 
eters to be determined (3 in our case). A good i ntroduction 
to this statistical approach can be found in iPress et al. (1992)). 
EWobs and EWsynth are respectively the EW of the observed 
spectrum and the EW of a synthetic spectrum, 5 red being com- 
puted for each synthetic spectrum of our 3D-grid. Lastly, cr cor- 
responds to the uncertainty on both EWobs and EWsymh such that 
cr^ = "^obs "'"'^synth- B^cause of the optimal normalization and the 
absence of noise in synthetic spectra, (Xsynth was neglected and 
we assumed that cr = iXobs, o-Qbs being described in Sect. 16.1.41 
and in Table|5]as crabs(EWobs)- The error bars on the parameters 
of the synthetic spectrum yielding the lowest value of S red were 
then investigated by means of statistical tools. The merit function 
5 red follows a distribution with v =N-k degrees of freedom 
if: i) we consider that the best model supplied by the minimum 
merit function is the right physical model representing the star 
and therefore that the corresponding stellar parameters are the 
true parameters of the star, ii) we assume that the perturbations 
on (EWobs - EWsynth) are independent and additive standard nor- 
mal random variables. Since the iXobs values are estimated from 
the local SNR around the selected lines and since the noise is 
poissonian and tends to a gaussian distribution in case of a great 
number of photons, which is the case here, the gaussian assump- 
tion is valid in this study. In this situation, the error bar corre- 
sponding to a choosen confidence interval x can be provided by 
the cumulative distribution function of the distribution with v 
degrees of freedom defined as: 



Dy(x) 



Jo I 



-dt - 



r(^) 



(5) 



where F and y denote the Gamma and lower incomplete Gamma 
function, respectively. D^ix) is equal to the probabilty P(5 < 
x) of the merit function S - 5 red x (N-k) to be within the 
confidence interval x. By misuse of langage, such confidence 
intervals on the S value are usually considered to be equivalent 
to confidence intervals on the resulting parameters. For the needs 
of our study, we did the same assumption. We were therefore 
searching for the interval A5 such that S + A5 < x, which was 
equivalent to search for A5red - A5/(N - k) such that 5red,min + 
AS' red ^ :c/(N - k). With V = 10 - 3 = 7 degrees of freedom we 
deduced that parameters of synthetic spectra included within an 
interval AS red of 1 . 1 7 , 2 . 04 and 3.12 from the minimum S red, min 
value were respectively 68.3%, 95.4% and 99.7% reliable. 

A preliminary computation of 5 red whithin the whole 3D- 
grid of synthetic spectra, whithout any constraint, yielded a 
Teff value of 8775 K, a logg value of 4.2 and a log Ape 
value of -4.25, i.e. [Fe/H]= H-0.29 (the solar abundance be- 
ing log(^^^ - -4.54). The corresponding error bars were 
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estimated using the method described above. For a confidence 
level of 68.3%, we found the following uncertainties: 525 K on 
Tetf, 0.7 on \ogg and 0.35 on log Ape. For a confidence level 
of 95.4%, we determined the uncertainties of 700 K on T^s and 
0.48 on log Ape, the small sensitivity to log g does not allow the 
determination of the 95.4% confidence level for \ogg within our 
parameter space. The results are shown in Fig.|7] 

What we pointed out at this stage were the limits of this first 
statistical computation due to the very low number of drastically 
selected lines used in the merit function minimization, that led to 
a significant Teff and log g degeneracy: numerous synthetic spec- 
tra yielded 5 red values within the chosen confidence intervals, 
and a much larger number of selected lines would have been 
necessary to converge reliably towards the model best represent- 
ing the observations with smaller error bars. This effect was also 
pointed out by the results of a jackknife resampling, which con- 
sisted in recomputing several times the merit function 5,ed, leav- 
ing out one after each other the equivalent width of one line. 
This statistical method revealed that the best parameters could 
change depending on the removed line within a AT^^ff of ;^2ook' 
a Alogg of and a Alog Ape of t'oog- This line-dependance 
of the result was attributed to the lack of selected lines too. To 
compensate for this lack, we confronted our preliminary result 
to the additional physical constraint brought by the excitation 
equilibrium: indeed, if T^ff is well determined, all spectral lines 
corresponding to a specified ion, and with a given excitation po- 
tential, should provide the same abundance determination. In a 
first step, for consistency, we only kept stellar model couples 
(Teff, \ogg) for which the equivalent width EW of all selected 
Fe I lines had an intersection with the curve of growth (EWsynth 
as a function of log Ape) within the range of our parameter space. 
In a second step we calculated for each selected model linear re- 
gression lines through the 8 data points logApe = f(x) (x being 
the excitation potential, usually expressed in eV) and kept only 
the stellar models yielding a close to zero slope regression line 
(in fact we calculated a standard deviation of the distribution of 
all slopes, and kept models with slopes 0. + 1 cr). Fig.[8]shows the 
fact that for each Teff a small range of log g and log Ape couples 
agreed with the excitation equilibrium within the error bar. 

Taking the excitation equilibrium into consideration, we ex- 
cluded models not respecting it and searched again for the min- 
imum value of S'red- We obtained: Ti^g- 8550 K, \ogg- 3.9 and 
log Ape - -4.38 . To determine the error bar on these parame- 
ters, we adopted the same approach as described previously in 
this section. It should be noted that log g is the less accurately 
determined stellar parameter because of the fainter sensibility 
of the selected lines to this parameter. The results can be seen 
in Fig. |9] and are summarized in Tabled The shape of the 2D- 
repartition od the 5 red values showed in Fig.|9]can be understood 
from Fig. |7] and Fig. [8] 



Table 6. Fundamental parameters of HD 104237 determined in 
this study. 



stellar 
par am. 


best model 
param. values 


68.3% 
conf. int. 


95.4% 
conf. int. 


(K ) 


8550 


+ 150 


±225 


logg 


3.9 


±0.3 


±0.4 


log Ape 


-4.38 


±0.19 


±0.26 


[Fe/H] 


+0.16 


±0.19 


±0.26 



In order to check the validity of our result, we plotted the 
abundances determined from the 8 Fe 1 selected lines as a func- 
tion of their excitation potential in the case of the model mini- 
mizing 5 red without any additional constraint and of the model 
minimizing 5 red with excitation equilibrium constraint. The re- 
sult can be seen in Fig.[TOl The linear regression across the points 
of the model with T^f^- 8775 K shows a negative slope, whereas 
the one across the points of the model with T^s- 8550 K shoes 
a flat regression line, which confirmed this model as the best so- 
lution. 



Excitation cquilibrijrn Vr\ 







1 1 1 

Model(Te-'f=8775K,logg= 1 1 -,=-l " 






* x _ 












Model(Te-'f=8550K,logg = 3.9,logAp,=-4.36j 



Exciia-jon potential (eVj 



Fig. 10. Excitation equilibrium of Fe I. In black: the model min- 
imizing 5 red with excitation equilibrium constraint. In grey: the 
model minimizing 5 red without any additional constraint. 



Another successful verification of the reliability of our re- 
sults consisted in verifying the iron ionization equilibrium, i.e. 
that an optimal model must give the same abundances derived 
from Fe I and Fe II. 




SOGO 



Fig. 11. Teff and log g pairs respecting the ionization equilibrium 
between Fe I and Fell lines in HD 104237. Stars: T^ff and logg 
values derived from the HD 104237 spectrum; the line shows a 
linear regression. 
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Fig. 7. 2D-repartition of the S red values as a function of Teff and log Ape ■ This graph has been built by considering successively every 
value of log g, from 3.5 to 4.5 by step of 0.1. The white X represents the best model, whose T^ff- 8775 K, logApe = -4.25 and 
log^= 4.2. The boundary 68.3%, 95.4% and 99.7% confidence intervals are respectively marked by the dotted lines and the black 
filled squares, the dashed lines and the grey filled circles and the black crosses. The grey cross represent the Tgif and log Ape ranges 
investigated in the present study from error box found in the literature, keeping in mind that log g values range from 3.5 to 4.5. The 
white filled square locates the best model with the additional excitation equihbrium constraint described in the continuation of the 
Sect.lO 



7. Discussion and Conclusions 

In this article we presented an in-depth analysis of the proto- 
type Herbig Ae star HD 104237. The star has been previously 
identified as a PMS 6 Scuti pulsator, and the comprehensive line 
profile analysis of the extended quasi-continuous high resolution 
spectroscopic time series obtained at SAAO from April 22'"^ - 
26"' 1999 has confirmed the main frequencies discovered by ra- 
dial velocity analysis in Bohm et al. (2004). For the first time, 
the less than 1.5% continuum level variations of the equiva- 
lent photospheric LSD profiles were studied and a direct con- 
firmation of the presence of at least one non radial pulsations 
could be performed. Based on a Fourier 2D analysis we iden- 
tified the dominant mode as an azimuthal order m = + 1, and 
a most likely { value of 1 or 2. Since HD 104237 is a moder- 
ate rotator with vsin/ of only 12+2 kms^'it will be difficult 
getting access to additional mode identifications in the future, a 
potential approach would be to organize multi-site continuous 
observations with significantly higher resolution echelle spec- 
trographs, ideally with R > 70000. On the mode identification 
aspect a modification of the FAMIAS code working with asym- 
metric photospheric profiles would be a great progress. 

In order to prepare a forthcoming asteroseismic modeling 
of this particular star we needed to determine very precisely 
its fundamental stellar parameters. Since HD 104237 is a mul- 
tiple system with a nearby spectroscopic binary companion, we 
decided to concentrate on the periastron night of 12''' of April 



2000, presen t in our data set of spectra acquired at SAAO (and 
described in iBohm et al.ll20d4l) . To combine all spectra of this 
particular night to one high SNR reference spectrum, we devel- 
oped an optimized spectral continuum determination tool. As a 
result, the reference spectrum had a SNR value of close to 400 
per resolved element at 550 nm. 

The detailed study of the fundamental stellar parameters has 
provided values of T^g = 8550 + 150 K, log^ = 3.9 + 0.3 and 
log Ape = 4.38 + 0.19 (i.e. [Fe/H]= H-0.16 + 0.19), the eiTor 
bars corresponding to the 68.3% confidence interval. A partic- 
ular effort was put on the statistically correct determination of 
the associated e rror bars . In a r eview of Teff and log g determi- 
nation methods. ISmallevI (l2005h mentioned optimal error bars to 
be of +100K for T^g, +0.2 for logg and of the order of 0.05 
to 0.1 in abundance. Due the high level of activity of HAeBe 
stars, and in the case our target star the spectroscopic binarity, 
it is very difficult to obtain a sufficiently large number of unpol- 
luted, exploitable photospheric spectral lines of a given ion; the 
achievable uncertainties remain therefore rather large. In addi- 
tion, a wrong estimate of the contribution of the secondary, both 
in spectral class and in luminosity, might lead to systematic er- 
rors which are very difficult to take into account. However, we 
estimate ffiese errors to be of second order in our study. 

The surface gravity being difficult to constrain in our spec- 
troscopic approach, the associated lumino sity ratio (base d on 
the empirical mass-luminosity relation by iMalkovl l2007h was 
determined to log (L^/L©) - 1.59 + 0.44, a value centered on 
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Fig. 8. (Teff, \ogg) pairs respecting the Fe I excitation equilibrium and associated log Ape ranges (dotted lines delimited by squares). 
Each color of square, from lightest grey to darkest grey, represents a value of log^ from 3.5 to 4.5. Grey crosses represent the 
subspace of models for which there exists an intersection between the EWobs line and the curve of growth for every one of the 8 
selected Fel lines. 



the result bv Ivan den A ncker et whose authors deter- 

mined log (L^/Lq) - 1.55^Q g^. On the contrary, our spectro- 
scopic effective temperature determination is significantly im- 
proved with Tpff = 8550 ±150 K, ag ain centered on the central 
value of Ivan den Ancker et al.l ( Il998l) . whose result was T^ff - 
8500+500 K. Our independent spectroscopic approach therefore 
tends t o confirm the previously pho tometrically determined val- 
ues bv lvan den Ancker et al.l (Il998h . Combining the two results 
would yield the following most likely fundamental parameters: 
log (U/Lo) = 1-55+°;°^, con-esponding to logg = 3.93 + 0.09, 
and Teff = 8550 + 150 K. Fig. [12] shows well the intersection 
betw een our new results and prev iously determined parameters 
from Ivan den Ancker et al.l (Il998h. At this stage it i s important 
to mention that error bars bv ^v an den Ancker et al.l (il998i) and 
iGradv et al.l (|2004) might be largely underestimated. 

lAcke & WaelkensI (|2004 determined an iron abundance of 
Fe ([Fe/H]=H-0.09 + 0.19. Our new result of the abundance, 
[Fe/H]= -hO.16 + 0.19, confirms t he conclusi o n that HP 104237 
might be slightly overabundant. IVick et all (1201 ll) calculated 
recently self-consistent stellar evolution models in this mass 
range including atomic diffusion; they concluded that for mass 
loss ranges above 10"'^ Moyr"' no surface abundance anoma- 
lies are expected. HD 104237 shows a very strong H„ line with 
a complex profile, which tends to indicate significantly highter 
mass loss rates, similar to other Herbig Ae stars showing typ- 
ical mass_Joss_jatesof_the order of 10"^ - lO^^Moyr ' (see 
e.g. iBouret & Catalal[l998h . Taking these new simulations into 
account, we conclude on a close to solar iron abundance for 
HD 104237, which is clearly within the error bars of both stud- 
ies. 



We want to mention at this stage that several assumptions 
have been made in this study, on the modeling side (LTE, plane- 
parrallel spectrum synthesis), on the contribution of the faint 
secondary (ratio of planck functions), on the threshold of line 
selection, but also the assumed photospheric line formation re- 
gion of our lines of interest, amongst others. Despite the fact 
that these residual errors should be negligible with respect to our 
study, further improvements on line selection and analysis, such 
as a combined treatment of blended lines, will lead to more and 
more constrained fundament al parameters in the future. 

A forthcoming article (iFumel et al.l 1201 If) will describe 
the asteroseismic modeling of the primary component of 
HD 104237. 
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Fig. 9. 2D-repartition of the Si-ad values as a function of Teff and logApe. The white X represents the best model, whose reff= 8550 K, 
logApe - -4.38 and \ogg= 3.9. Confidence interval symbols are the same as those used in Fig|7] Grey crosses have the same 
meaning as in Fig. [8] 
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